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I. INTRODUCTION 

The power counting renormalizable non-relativistic theory of gravity recently proposed by Hof ava [l| , is a theory 
of gravity in which general covariance is partially abandoned in favor of renormalizability. A state of the theory is 
defined by a four-dimensional manifold A4 equipped with a three dimensional foliation J 7 , with a pseudo-Riemannian 
structure defined by an Euclidean three dimensional metric in each slice of the foliation gij (x, t) , a shift vector N 1 (x, t) 
and a lapse function N(x,t). This structure can be encoded in the ADM-decomposed metric 

ds 2 = -N 2 (x, t)dt 2 + gij{x, t) [dx 1 + N\S, t)dt) (dx j + N j (x*, i)dt) . (1) 

The dynamics for the set (M , J-, g^ , N t , N) is defined as being gauge invariant with respect to foliation-preserving 
diffcomorfisms, and having a UV fixed point at z = 3, where the dynamical critical exponent z is defined as the 
scaling dimension of time as compared to that of space directions [x] = —1, [t] = — z. This choice leads to power 
counting renormalizability of the theory in the UV. To the resulting action one may add relevant deformations given 
by operators of lower dimensions, that lead the theory to a IR fixed point with z = 1, in which symmetry between 
space and time is restored, and thus a generally covariant theory may emerge. 

In order to have control on the number of terms arising as possible potential terms, one may impose the so-called 
detailed balance condition: the potential term in the action for 3+1 dimensional non-relativistic gravity is built 
from the square of the functional derivative of a suitable action for Euclidean three-dimensional gravity (here three- 
dimensional indices are contracted with the inverse De Witt metric). Condensed matter experience on this kind of 
construction tells us that the higher dimensional theory satisfying the detailed balance condition inherits the quantum 
properties of the lower dimensional one. It has to be noted that the theory is still well defined even when detailed 
balance condition is broken softly, in the sense of adding relevant operators of dimension lower than that of the 
operators appearing at the short distance fixed point z = 3. With such a deformation, in the UV the theory satisfies 
detailed balance, while in the IR the theory flows to a z = 1 fixed point. 

We won't go through the above described steps in more detail, but state the resulting action that will be relevant 
to our purposes. The interested reader can refer to the original paper [![. The action for non-relativistic gravity 
satisfying the detailed balance condition can be written as 

( * (Kl] K- \K 2 ) + *V(jWfi - 3A 2 w -) + *V(l- 4A) fi2 _ # f ^ _ f &j _ f^ Rij \ \ >(JJ) 



8(1 -3A) 32(1 -3A) 2w 4 \ J 2 J J\ 2 

Here A\y, k, A, fj, and w are arbitrary couplings, R 7 Rij,Cij and K%j are the scalar curvature, the Ricci tensor, the 
Cotton- York tensor and the extrinsic curvature respectively of the three-dimensional sections of the foliation. The 
dynamics in the infrared is controlled by the first two terms, and then, if A = 1, general relativity is recovered. On 
the other hand, in the UV the third and fourth terms become dominant, and the anisotropy between space and time 
is explicit. 

This action is invariant under foliation preserving diffeomorphisms, namely under changes of coordinates of the 
form 

x*' =x i \x j ,t) 

t'=t'(t) (3) 
under which the spatial metric, shift vector and lapse function transform as 

g' ij (x i ',t') = d i ,x k d j ,x l g kl (x r ,t) 
N i '(x i ',t') = d k x i ' d t >tN k (x j ,t) 

N'tf' \t') = d t dN{x\t) (4) 

Eqs. ([3]) and (|4]) ensure that, if the lapse function is initially chosen to be independent of the space coordinates in a 
given coordinate system, it cannot be turned into a space-dependent form by a change of coordinates. In other words, 
the space-independence of the lapse function is a covariant statement. This implies the existence of two possible 
versions of Hofava gravity, a "projectable theory" in which the lapse function is space-independent N — N(t), and a 
"non-projectable theory" in which the lapse function is allowed to depend on space N = N(x 1 ,t). 
In the non-projectable case, the equations of motion obtained by varying the above action are 

-(K-K** XK 2 ) ^Ml^) *V(1-4A) 2 k* „ 

k2 {K v K 8(1 -3A) " 32(1 - 3A) R + 2w± ^ ~°' (5) 

V k {K u -\Kg u ) = Q, (6) 
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where 



and 



A _ . K V A 1V p(3) , « 2 M 2 (1 - 4A) (4 ) _ ( 5 ) _ jf_ (6) ,„s 

K 2 « K 2 ij " l "8(l-3A) y 32(1 -3A) « 4w 2 « 2w 4 « kU 



= Cy - — — -Rij , (8) 



£# } = NiS7 k K k j + NjV k K k i - K k iS7jN k - K^Nk - N k V k K K 
—2NK ik Kj k - X -NK kl K kl 9ij + NKKy + K %3 , 



i j 



\nK 2 9ij + NidjK + NjdiK - N k {d k K) 9ij + K 9ij , 



E\f = N(Rij - \R m + \k w9ij ) - (ViVj - 9 ijV fe V fc )iV . 
1 



E\f = NR{2R ij - -R 9lj ) - 2(ViVj - 9l] V k V k )(NR) , 



E if = Vk[Vi(NZ k i) + V l {NZ k J )} - V k V k {NZ tj ) - V k V,{NZ M ) m , 



E\f = - X -NZ kf Z kl 9l3 + 2NZ ik Z k ~ N{Z lk C k + Z jk C k ) + NZ u C H 9ij 



— ■^7 k[Ne mk (Z mi Rj£ + Z m jRig)] + -R T ' \ V n [Ne mke (Z mi 
-^7 n [NZ r :e mkl { 9kl R ol + 9k3 R u )\ - ^ n W n W k [Ne mk "(Z mi9]e + Z mj9u )] 
+ \v n [V l V k (NZ^e mke )g 3e + V k {N Z^ H ) 9u ] 

+lv e [ViV k {NZ mj e mkl ) + V 3 V k {NZ ml e mM )} - V ' n V e V ' k (N Z"e mki ) 9ij . (9) 

In the projectable case, eq.© is replaced by its spatial integral. 

Since the original proposal of [l[ , there has been a growing number of research papers in the area. Formal devel- 
opments were presented in some spherically symmetric solutions were presented in [Io|-[l6j. rotating solutions 
were studied in [13], string- like ansatze were investigated in [l8j]- [l9j], toroidal solutions were found in (20|, gravita- 
tional waves were studied in [2l| and (22[, cosmological implications were investigated in [23l]-[35|h and interesting 
features of field theory in curved space and black hole physics were presented in (36|-[42j]. In |43j|-[46j potentially 
harmful instabilities were pointed out, originated in the additional scalar graviton mode, that propagates in virtue 
of the reduced gauge symmetry. Moreover, there it was shown that the extra mode become strongly coupled at the 
infrared in nontrivial backgrounds. To cure these problems, a so called "healty extension" of the non-pro jectable 
theory was proposed, in which additional terms containing derivatives of the lapse function were included in the ac- 
tion, which has the effect of eliminating the instabilities [471 ] . Alternatively, a covariant theory whose partially gauge 
fixed version reproduces non-projectable Hof ava dynamics was developed, in which the extra mode was shown to be 
harmless [48[ . In [49| , a reinterpretation of a secondary constraint that appears when A ^ 1 in the infrared limit of the 
non-projectable theory, leads to a new counting of degrees of freedom in which the extra mode is not present. Finally, 
in |50| | an additional U(l) gauge symmetry was introduced, that kills the scalar graviton avoiding the aforementioned 
problems. 



In Einstein gravity, it is easy to prove that no non-trivial solution with the symmetry of a domain wall exist in the 
absence of matter. Indeed, the only solution of the equations of motion compatible with a smooth and flat domain 
wall ansatz is that of an AdS/Mikowski spacetime, depending on the cosmological constant. In Hof ava gravity on the 
other hand, the terms containing higher spatial derivatives could in principle play the role of a matter contribution, 
allowing for the existence of non-trivial domain wall solutions in vacuum. This is one of the motivations of our work. 

In Einstein gravity, the knowledge of the cosmological solution corresponding to a given kind of matter can be used 
to obtain a domain wall solution through the so-called domain wall/cosmology correspondence (5l| . Indeed, given a 
metric with the Friedmann-Lemaitre-Robertson- Walker form, it can be mapped to a domain wall ansatz via suitably 
defined Wick rotations of the coordinates. Moreover, such transformation maps the Friedman equations into the 
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equation of motion corresponding to the domain wall, ensuring that cosmological solutions are mapped into domain 
wall solutions. In Hofava gravity on the other hand, the situation is very different. The essential anisotropy between 
space and time present in the theory is an obstacle for the domain wall/cosmology correspondence to work. First, a 
cosmological ansatz has a single independent function that can be identified with the scale factor while, as we will see 
bellow in further detail, a domain wall ansatz has in principle two independent functions. Moreover, the anisotropy 
between space and time implies that the equations of motion for a domain wall are not mapped under Wick rotation 
into the Friedmann-like equations for a cosmological ansatz. In consequence, the large amount of research regarding 
cosmological solutions of Hofava gravity [23|-|35| gives no clue on the form of the domain wall solutions of the theory, 
making the study of such solutions a subject of independent inquiry. This is a second motivation for the present 
paper. 

In this paper we start the investigation of domain wall solutions of Hofava theory. The simplest possible setup 
being that of a purely gravitational theory, we will not include any matter degree of freedom in our equations. In 
attention to the the fact that all the solutions to the aforementioned controversy about the extra scalar mode that were 
proposed in the literature correspond to modifications of the non-projectable theory, we will limit our investigation 
to that case. Since in the absence of matter, there is no possible Z2 symmetry, to be broken differently at each side of 
the wall forming "domains", we will sometimes use the somewhat more accurate name "membrane" for our solutions. 

II. DOMAIN WALL SOLUTIONS 

This paper deals with the issue of domain wall solutions of Hofava theory. In the present non-relativistic context, a 
flat domain wall solution is defined as a solution having translational and rotational symmetry in two dimensions, i.e. 
being invariant with respect to the ISO(2) x K group of transformations. The symmetry of the solution is reduced 
with respect to that of a relativistic domain wall ISO(2, 1) because of the non-relativistic nature of Hofava theory. 
An ansatz that preserves such symmetry can be easily written as 

ds 2 = -e v ^dt 2 + e U ^ (dx 2 + dy 2 ) + dz 2 . (10) 

Note that, in virtue of the reduced symmetry of the theory, there is no set of coordinates in which U(z) = V(z) 
as it would happen in Einstein theory, and then the ansatz has two independent functions to be determined by the 
equations of motion. The variables x, y can be chosen as describing a two-torus T2 of volume V xy or a two- plane R 2 
(that can be considered as the infinite V xy limit). In what follow we will call "lapse function" to e v ^ and "spatial 
volume function" to e u ^ z K 

Replacing the ansatz in the above equations (JSJ) - (O we get the following two independent equations of motion 

(4A^ + U' 2 ) (3 (4A W + U' 2 ) + 8U") - 8(A - 1)U" 2 = , (11) 

48A 2 y + ((4V' - U') U' 2 + 8A W {V + 2V')) U' + 8(A - 1) ({U" - 2U' 2 - U'V) U" - 2U'U {3) ^ = . (12) 

Note that in the above equations the constants /z, k and w do not appear. This can be traced back to the action 
@ or to the equations of motion (J5j- in which, when the extrinsic curvature and the Cotton- York tensor CV,- 
vanish, as it happens for our ansatz, the w constant cancels and the product K 2 fi 2 can be factored out. As can be seen 
in the second equation, the case A = 1 is special in the fact that the total differential order of the system is reduced. 
As we will see, this property manifest itself in a non-analyticity of the solutions as functions of the parameter A. 

A. Solutions with Aw < 

1. Solutions with A = 1 

We first fix our attention in the case A = 1. Solving cq. tfTTj) for U(z) and replacing the solution into eq.(Q~2|) to get 
V(z), we obtain the corresponding solution. It reads 

e v.M = cosh (hE^ {z Zo) y\ mh (hE*Z. {z - Zo ij , (13) 
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where z a is a constant of integration, another two constants of integration have been reabsorbed in the definition of t 
and of x, y. We see that the solution is Z 2 symmetric, and that its asymptotic form for z — > ±00 is that of an AdS 
spacetime 



In order to have a physical interpretation of this solution, we evaluate some of the observable scalars of the theory. 
We start with the spacetime curvature, that reads 



# = ^fll + 5taai^V^(z-J ) — ► 12k w , (15) 
4 \ \4 J J z->±oo 

We see that it approaches a constant value when z — » ±00, as may have been expected from is asymptotic AdS 
form. More interestingly, it shows a peak at z = z a . A low energy observer, armed only with the tools of Einstein 
gravity, would conclude that some kind of matter with a positive energy density has to be sitting there, to partially 
cancel the contribution of the cosmological constant. By observing that the total spatial volume of the z — z a slice is 
proportional to V xy , he/she would identify the matter distribution as a "membrane" located at z — z Q . Nevertheless, 
from our privileged high energy point of view, we know that what we have is a purely gravitational soliton, since 
no additional matter has been added to Hofava theory. With this at hand, we will refer to our solution as a purely 
gravitational membrane located at z — z Q . 

Since Hofava theory distinguishes explicitly time from space, a separated invariant of interest is the space curvature, 
that in our case is given by 



i?( 3 ) = 3A W 1^1 + tanh [^V^A w(z - z )\ — > 6A W . (16) 

Again it presents a peak at z — z Ql that reinforces the identification of that point as the location of the membrane. 
Plots of the space and spacetime curvatures have been included in FigJT] 

This is all we can say about our solution from a high energy point of view. A low energy observer on the other 
hand, will experience a world in which the velocity of light constitutes an upper bound to the propagation of signals, 
and where in consequence the light-cones of the metric determine the causal structure. From that perspective he/she 
may wonder about the presence of horizons and/or hidden regions. To investigate such possibility, we calculate the 
time taken by a light signal emitted for an observer that sits at an arbitrary point z\ to reach a different point z 2 . It 

is given by At(z2,zi) = e (£(22) — K z i))i where 



t(z) — dz e 2 



2 1 

3' 3 



(17) 



3V^A^ 

where we have defined the notation £(z) = sinh 2 (3\/ Aw (z ~ z Q )/4) , and B x [p, q] is the Euler incomplete Beta function. 

When z — > z Q , we have £(z) — > and B_i/£( z )[2/3, 1/3] ~ — (— l)s log[£(z)] diverges. This implies that the time taken 
by a light signal to reach the slice at the center of the geometry z — z Q is infinite, which allows us to identify such 
slice as a horizon. 

At this point a remark is in order: following the intuition gained on the study covariant theories, one may think 
that a change of variables such like r 2 = e Uo ^ can be used to analytically extend the solution beyond the horizon, to 
the region r < 1. In that region, the t variable would become spacelike, while the r variable would become timelike, 
implying that the equations of motion satisfied by the metric (fTT|) - ([T2)) contain third order time derivatives. But, by 
construction, Hofava action in eq.@ cannot give rise to equations of motion with more than two time derivatives. In 
consequence, in the region r < 1 the analytically continued metric is not a solution of action @. This implies that 
such an analytic continuation is not allowed in Hofava theory. 

To clarify the causal structure of a metric, the usual technique in relativistic theories is to isolate two of the 
dimensions of the metric (time and one of the space dimensions, in our case it would be z), and then perform a change 
of coordinates in order to bring the metric of the z,t slices into a conformally flat form, with the new coordinates 
running in a compact interval. The resulting Penrose diagram is then easily drawn, by stripping the conformal factor. 
Indeed, if one defines new coordinates as 

u± = Sgn(z) Tanh (e~ t(2)±t ) (18) 




z — 




FIG. 1. Plots of the lapse function goo = —e"^ z ' (up-left) and the spatial volume function g xx = g yv — e u °^ z ' (up-right), 
spacetime curvature R^'(z) (down-left) and space curvature R^ 3 '(z) (down-right), for Aw < and A = 1 



the resulting metric reads 



where now the u± variables are defined inside the square [—1,1] x [—1,1], in which light rays propagate along the u± 
directions. The the conformal factor of the u+,U- slice can now be stripped, which allows for the construction of the 
corresponding Penrose diagram. It is given in Figf?l 

Nevertheless, at this point a very important remark is in order: in Hofava theory, changes of coordinates that mix 
time with space do not leave the action invariant, i.e. they are not symmetries of the theory. In consequence metric 
(|19|). obtained from (fT0|) by the illegal change of coordinates (fT8|) . is not guaranteed to be a solution of the equations 
of motion. For that reason the Penrose diagram in Figf?] can be considered only as an approximation valid at low 
energies, where general relativity is recovered. In order to have a description of the causal structure valid at any 
energy, the most one can do is to draw a diagram showing the light cones of the metric in the original coordinates, or 
in some new coordinates obtained from them via foliation preserving diffeomorphisms. We included such diagram in 



ds 2 



e 2t(z)+V{z) 



^(dx 2 +dy 2 ) 



(19) 



du+du- + e 



FigH 
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2. Solutions with A / 1 

Let us next assume that A ^ 1. Then solving eq. (fTTj) for U{z) and replacing the solution into ea. (lT2l to get V(z) 
we have 

2p±(A) 



\ p±W J 



e v *W = (^cosh (^(* - *)) W±W "" -nh (^(* " *■) ) ) = N* , (20) 
with 

2(A-1) , N 

and where again z D is a constant of integration and another two constants have been reabsorbed in the definition of 
t and of x, y. Here the subindex ± indicates that we have two different branches of solutions according to the choice 
of sign in p±(X). Note that for A < 1/3 no real solution exists, so in the remaining of this subsection we will focus on 
the region of parameters 1 7^ A > 1/3. 

Both ± solutions are Z2 symmetric and centered at z . To explore the asymptotic behavior we take z — > ±00 to 
have 



e U±{z) x e 2V-A w sign(p±(A)) \e-z„\ 



gV±(z) ^ e 2z ± (A) v / -Aiv sign(p±(A)) \z-z a \ _ (22) 

where 

4(2 P± (A)-A-1) 
P±(A)(3p±(A) - 2) 

Here we see that our metric corresponds to an asymptotically Lifshitz spacetime, similar to those studied in [52j . 
whose scaling exponent is given by z±(A). At this point, it is convenient to stress that such scaling exponent is not 
in principle related with the dynamical critical exponent z = 3 of Hofava theory. 

In the solution with the + sign, p+(A) is always positive, and then the spatial volume function is exponentially 
growing when z — > ±00 independently of the value of the parameter A. The lapse function instead grows exponentially 
whenever z+(A) > (A < 3) and decreases exponentially for z+(A) < (A > 3). In the limiting case z + (A) = 
(A = 3) the function asymptotes a constant value. Plots of the corresponding solutions for different values of A can 
be seen in Figf2J 

On the other hand, in the solution with the — sign, the spatial volume function is exponentially decreasing for 
z —> ±00 when p_(A) < (A > 1), and exponentially growing for p_(A) > (A < 1). Additionally, z_(A) is always 
positive, which implies that the lapse function grows exponentially when z — > ±00 for any value of A. Plots of the 
corresponding solutions for different values of A can be seen in FigJ3] 

Again, to have a physical interpretation, we evaluate the spacetime curvature scalar for this solutions. It reads 

m = i?^(A) + (i£>(A) - R^W) sech ( VESip^ 2 . (24) 

Its asymptotic value is consistent with a Lifshitz spacetime and it is given by 

(4) _ -2A W (-16(1 + A) 2 + p±(A)(48(l + A) - P± (A)(20 - 24A + 3 P± (A)(4 + 9 P± (A))))) 

(2 - ■ip±{X)) 2 p±{Xy 

As before, a peak at z — z a appears, with the value 

R(4)m R ( 4)m 2A W (8A(1 + 2A) + P± (A)(-4(1 + 3A) + P± (A)(14 - 24A + 3p±(A)(-8 + 9p±(A))))) 

R * (A) = RJ (A) (2-3 P± (A))2p ± (A)2 ' (26) 
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FIG. 2. Plots of the lapse function goo = —e v +^ z ' (up- left) and the spatial volume function g xx — g yy — e u +^ (up-right), 
spacetime curvature i? (4) (z) (down-left) and space curvature R!- 3 \z) (down-right) for Aw < 0, for different values of A. Notice 
the change of the asymptotic behavior of e v+( - z ' at A = 3 in the upper left plot. 



Such peak, together with the observation that the spatial volume of the z — z a slice is proportional to V xy , allows 
us to make the same interpretation as before, referring to our solution as a purely gravitational membrane sitting at 
z = z . This interpretation is reinforced by evaluating the space curvature scalar, that is given by 



flffl =2A W U+ sech ( ^f- Zo) ) -^ = 6A W , (27) 

and that also shows a peak at z = z . As can be easily seen in the above equation, the space curvature at infinity is 
independent of A. On the other hand its sign at z a depends on A, for the solution with the + sign it is always positive, 
while for the solution with the — sign it is positive for A > 1 and negative otherwise. 

Plots of both spacetime and space curvatures for each kind of solution ± were included in Figs. [5] and [31 

To investigate the presence of a horizon at z = z where the lapse function vanishes, again we evaluate the time 

taken by a light signal emitted for an observer that sits at an arbitrary point z = Z\ to reach a different point z — z 2 ■ 

y±(*i) 

It is now given by Ar(z2,zi) = e 2 (£(2:2) — ^( z i)): with 

f z v ±w 
t(z) = dz e 2 = 

= 1 ] - X P±( > B_ imz) [l+q ± (X),- q± (X)} (28) 
V -Aw 

where we used a generalization of our previous shorthand notation £(z) = sinh 2 (y/ k\v(z — z )/p±(\)}. The number 
q±{\) has been defined as 



(29) 






FIG. 3. Plots of the lapse function goo = 



(up-left) and th spatial volume function 



— 9yy 



(up-right), 



spacetime curvature R (4) (z) (down-left) and space curvature R^(z) (down-right) for Aw < for different values of A. Notice 
the change of the asymptotic behavior of e u -^ at A = 1 in the upper right plot. 



By choosing z = z D , we have £(z) = and the Euler incomplete Beta function B_i/£( z )[l + <?±(A), — <?±(A)] ~ 
(— 1)9±( A ) log[£(;z)] diverges. This implies that At(z ,zi) is infinite, which allows us to identify the center of the 
geometry z = z as a horizon. The plot of the light cones corresponding to the present solution is given at Fig. 21 
where its approximate Penrose diagram was also included. 

Notice that these solutions are analytic in the parameter A in all its range except at A = 1. There, the solution 
with the + sign is analytic and approaches the solution for A = 1 given in (|13[) . On the other hand, the solution with 
the — sign is non-analytic at that point, since p_ (A) vanishes there. 



3. Degenerate solutions 

In the present Aw < case, eas.(fTT ]) -([T2 |) allow for an additional (degenerate) branch of solutions for any value of 
A. It is given by 



e U d {z) _ e 2y^W (z-Zo) 



,V d (z) _ 



arbitrary function . 



(30) 



The existence of such an infinite branch of solutions, for which the lapse function is not determined, was alr eady 
pointed out in the spherically symmetric case in [Io| and in the warped BTZ string context in 18]. As shown in 10], 
this degeneracy is closely related to the detailed balance condition, and is lifted by an small violation of it. 
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a 



Horizon 




FIG. 4. Left: plot of the light cones corresponding to the solutions with A < 0. Right: Penrose diagram of the same solution. 
Region I (II) corresponds to z > (z < 0), while regions III and IV are absent on the original solution. Since the changes 
of coordinates needed to strip the conformal factor mix time with space, the Penrose diagram must be considered only as an 
approximate description, valid only at low energies. 

B. Solutions with Aw > 

1. Solutions with A = 1 

Proceding as in the previous sections, we first assume that A = 1. The corresponding solution reads 



e u 'W = cos 



e v ° (z) = cos 




(z - z Q ) 



8/3 



{z-z )\ sin ( ^^f^ -{z - z ) 



(31) 



where, as before, z is a constant of integration, the solution being Z2 symmetric around z = z a . The lapse function 
is real in the region \z — z a \ < 2tt/3\/A w and imaginary otherwise. Then the surfaces z — z ± 2n/2>^/X w define the 
boundaries of spacetime. At those boundaries, the lapse function and the function e Uo ^ vanish. 

To have a physical interpretation of the solution, we evaluate the curvature scalar, that takes the form 



R {4) = ^A w I 11 - 5tan (^y/X^(z 




-co . 



Z— ¥Z a ±- 



(32) 



It takes a finite positive value proportional to Aw at z — z a , and it blows up to negative values at the boundaries. 
From the point of view of a low energy observer, that interprets the solution in the light of Einstein theory, some 
kind of negative energy density has to be localized close to each boundary, in order to cancel the contribution of the 
positive cosmological constant, resulting into a negative spacetime curvature. From our high energy point of view on 
the other hand, we know that we are in the presence of a purely gravitational soliton. Since e Uo ^ vanishes at the 
boundaries, the spatial volume of the slices z — z ± 2-k ji\fA w is zero, implying that the singularities are point-like. 
On the other hand, all the remaining slices z 7^ z a ± 27r/3\/A\v have a finite non-zero volume proportional to V xy . 
Then the interpretation of this solution is that of two point like singularities sitting at the poles of the geometry 
z = z a ± 2tt/3VAw and separated by a finite physical distance 4ir/3\/A w . The space curvature confirm such view 



-3 A 



w 



-1 + tan ( -^J A w (z - z ) 



-00 . 



Z— >Zo±~ 



(33) 



being again singular at the boundaries, while at the center it is positive and finite, proportional to Aw. 



_ e v+W 



z — z a 



e V+{z) 
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FIG. 5. Plots of the lapse function <?oo = — e (up-left) and the spatial volume function <7 ra = = e c/ °' z ' (up-right), 
spacetime curvature (z) (down-left) and space curvature -R' 3 '(z) (down-right) for Aw > and A = 1 



Plots of the solutions together with their spatial and spacetime curvatures can be seen in FigfSJ 



To investigate whether the singularity we just found is hidden behind a horizon, we proceed as in the previous 

V (zi) 

section. The time taken by a light signal emitted at z = z\ to reach z = Z2 is Ar(z2, z i) = e s (t( z 2) — t( z i)) where 
now 



t{z) 



dz e 2 = 

1 (- /r _ 1 (2 X 1/3 (z) + l 



2V3cot~ 



V3 



In 1 



3 X V3 (JS ) 



X 2 / 3 (z)+x 1/3 (z) + l 



(34) 



where we used the notation x( z ) — cos 2 (3y/Aw{z — z )/i). By choosing z = z , we have x( z ) = 1 an d the logarithm 
diverges. This implies that the slice at the center of the geometry z — z a is a horizon. On the other hand, by choosing 
z = z a ± 27r/3\/Avv we have x( z ) = Oj an( i the integral is finite. This implies that there are no horizons sitting at 
the singularities. In conclusion, from the point of view of any observer, one of the singularities is beyond a horizon 
and the other is naked. We included the corresponding plot of the light cones in FigJ5] and the approximate Penrose 
diagram in FigJHl 
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2. Solutions with A / 1 



Next, we move to the case A ^ 1 . Then, solving eq. flTTt for U{z) and replacing the solution into eq. (|T2"j) to get 
V(z), we have 



e^±W = cos (z - z a ) 



P±(A) 



(2 - Zo) 



2p±(A) 



5p_^(A)-4A-2 
3p ± (A)-2 



sin I — — (2 - z ) 



p±W 



(35) 



where, as before, z a is a constant of integration, and p±(X) is given by expression (|21[) . Again the solution does not 
exists for A < 1/3, while for A > 1/3 we get two branches of solutions identified with the ± subindices according to 
the choice of sign in the p±(\). 



v+M 



-e v + 



U+(z) 



e u + 




R {i \z) RW(z) 




FIG. 6. Plots of the lapse function goo = — e v+l - z ' (up-left) and the spatial volume g xx = g yy = e u +^ z ' ) (up-right), spacetime 
curvature R w (z) (down-left) and space curvature R^(z) (down-right) for Aw > 0, for different values of A. Notice the change 
of behavior of the lapse function at the horizons at A = 3, in the top left figure. 

For both branches, the lapse function is real in the region \z — z Q \ < irp±(\)/2\/ k\y and imaginary otherwise. Then 
the surfaces z = z a ± irp±{\) /2V define the boundaries of space. At the boundary, the behavior of the solution 
depends on the value of A. For the solution with the + sign, the exponent of the cosine in the lapse function changes 
sign at A = 3, implying that for A < 3 the lapse function vanishes at the boundary, while for A > 3 it diverges. On the 
other hand, for the solution with the minus sign, the exponent of the cosine in the spatial volume function changes 
sign when A = 1, resulting in a solution whose spatial volume vanishes at the boundary for A < 1 and diverges for 
A > 1. 

The spacetime curvature for the above solutions read 

RW = flW(A) + (4 4) (A) - £ (4) (A)) sec ( ~ Z ° ] V , (36) 

V p±W J 
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FIG. 7. Plots of the lapse function goo = —e v -^ z ' (up-left) and the spatial volume function g xx = g vy — e u ~^ z ' (up-right), 
spacetime curvature R^ A \z) (down- left) and space curvature R^ 3 \z) (down-right) for Aw > 0, for different values of A. Notice 
the change of behavior of the lapse function at the horizons at A = 1, in the top right figure. 



and 



where 

- (4) -2A W (-16(1 + A) 2 + P± (A)(48(1 + A) - P± (A)(20 - 24A + 3 P± (A)(4 + 9 P± (A))))) 

H " A J = To — q — TTxvi — 7Tv2 ' \ 67 ) 

R( 4) (A) = ^(4) (A) _ 2A W (8A(1 + 2A) + P± (A)(-4(1 + 3A) + P± (A)(14 - 24A + 3 P± (A)(-8 + 9 P± (A))))) 

(2 - 3p±(A)) 2 p±(A) 2 

Again it is clear that the solution is singular at the boundaries, where the function sec (y/— A\y(z — z )/p±(A)) 
diverges. These singularities are separated by a finite distance irp±(A) / \J Aw ■ 

Since the spatial volume function vanishes at the boundaries whenever p± (A) is positive, in such case the singularities 
will be point-like. This happens for any A in the + branch of solutions, and for A < 1 in the — branch. In the rest 
of the — branch, that is for A > 1, the volume of the boundaries diverges, implying that the singularities can be 
interpreted as two parallel purely gravitational membranes, sitting at the boundaries of spacetime. 

Plots of the solutions with the + sign with the corresponding spacetime and space curvatures for different values of 
A can be seen in FigJSJ On the other hand, plots of the solutions with the — sign with the corresponding spacetime 
and space curvatures for different values of A can be seen in Fig(7J 

The time taken by a light signal emitted from an observer sitting at z = z\ to reach a point z — z-i is given by 

At(z 2 , zi) = e a (t(z 2 ) - t(zi)) with 

f z v ± ( *> 
t(z) = dze 2 = 

'" x q±W (z)2F 1 {l, q± (\);l + q ± (\); X (z)} (39) 
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where the notation is now x( z ) = cos 2 (y/Aw(z — z )/p±{\)), the symbol 2-F1 representing the hypergeometric func- 
tion. Since such function diverges at xi z ) = 1> that is at z — z a , the time taken by a light signal to reach such slice 
is infinite. In consequence, as in the case A = 1, there will be a horizon at the center of space. On the other hand, 
the prefactor x q± ^H z ) vanishes at z — z a ± 7rp±(A) /2V A\y whenever the exponent q±(X) is positive, and diverges 
whenever it is negative. This implies that, as in the case A = 1, there will be no horizons at the singularities for 
A > 1/2 in the + branch. For the rest of the + branch, that is for A < 1/2, as well as for all values of A in the — 
branch, there will be a horizon sitting at each of the singularities. An observer at any point of the bulk will regard 
one or both of the singularities as hidden behind horizons according to the value of A and to the branch he/she is 
probing. 

As conclusion, the solutions of the + branch will consist in two point-like singularities separated by a finite distance, 
which arc hidden behind horizons whenever A < 1/2. On the other hand, the solutions of the — branch will consist 
for A < 1 in two point like singularities, and for A > 1 in two parallel membranes, in both cases the singularities being 
hidden behind horizons. Plots of the light cones in both cases can be seen in Fig|8l the approximate Penrose diagrams 
can be seen in FigJ5J 




FIG. 8. Plots of the light cones for the Aw > solutions. On the left we see the plot corresponding to the case in which there 
is a single horizon at the center of space, that is A > 1/2 in the + branch, which includes the nondegenerate solution with 
A = 1. On the right, the case in which the singularities are hidden behind horizons, that is A < 1/2 in the + branch and all 
values of A in the — branch. 



III. DISCUSSION 



We explored membrane solutions to Hof ava nonrelativistic theory of gravity when the detailed balance condition 
is satisfied. We found that for arbitrary values of the parameter A > 1/3, branches of membrane solutions exist. 

For the particular case A = 1 there is a single branch that corresponds to a Z2 symmetric spacetime. In the case 
of negative cosmological constant, the spacetime and spatial curvatures have a peak at z = z Q , that allows us to 
identify the solution as a membrane sitting at the center of space. There is a horizon at the membrane location. The 
curvatures become proportional to An/ at large distances, where the metric correspond to an AdS spacetime. On the 
other hand, in the case of positive cosmological constant, the space has boundaries at \z — z Q \ — 2n/3y/A w - beyond 
which the metric becomes complex. The spacetime and space curvatures diverge at the boundaries. The spatial 
volume of the slices containing the singularities vanishes, which implies that the singularities are point like. There 
is a horizon at the center of space, implying that from the point of view of any observer, one of the singularities is 
hidden beyond a horizon an the other is naked. On the other hand, both curvatures are finite and positive at the 
intermediate space. 

For generic values of A in the region 1/3 < A 7^ 1, two branches appear that correspond to I2 symmetric solutions. 
In the case of negative cosmological constant, the space is unbounded and the space curvature asymptotes a constant 
value proportional to Ayy. The same is true for the spacetime curvature, but with a proportionality factor that is 
a function of A. The asymptotic metric corresponds to a Lifshitz spacetime. Again, the curvatures being finitely 
peaked at the center of space, we interpret the solution as representing a membrane sitting there. The curvatures 
at the center depend on the value of A. There is a horizon at the membrane. On the other hand, in the case of 
positive cosmological constant, the space is bounded and the curvatures are singular at the boundaries. The spatial 
volume of those boundaries vanishes for A < 1 in the — branch and for all values of A in the + branch, implying that 
the singularities are point like. On the other hand for A > 1 in the — branch the area of the slices containing the 
singularities diverges, allowing the interpretation of them as membranes. As in the A = 1 case, there is a horizon 
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FIG. 9. Penrose diagrams of the A > solutions. The left (right) figure correspond to A > 1/2 in the + branch (A < 1/2 in 
the + branch and all values of A in the — branch). Regions I (II) correspond to z > (z < 0), while regions III and IV are 
absent in the original solution. In the left plot the singularities are naked, while in the right they are hidden behind horizons. 
Again, it must be kept in mind that these diagrams are only valid as a low energy approximation, since in constructing them 
we used changes of coordinates that mix time with space. 

at the center of space for any value of A in both branches. For A > 1/2 in + branch, there is no horizon at the 
boundaries, while for A < 1/2 in the + branch and for all values of A in the — branch the singularities are hidden 
behind horizons. 

It should be kept in mind that our solutions are purely gravitational solitons, since no additional matter terms 
have been added to Hofava action. Nevertheless, our nomenclature was inspired in the point of view of a low energy 
observer, according to whom the dynamics of gravity is totally covariant and given by Einstein theory. He/she would 
necessarily interpret our gravitational domain wall as originated on some kind of membrane-like matter sources. 

For any value of A, and additional branch exists. It is degenerate, in the sense that the lapse function is completely 
undetermined by the equations of motion. This behavior has been reported before in the case of spherically symmetric 
[To| and warped BTZ string [l8[ solutions, and was related with the detailed balance condition [10l |. 

It is interesting to note that the solutions behave analytically in the parameter A for all values of A ^ 1. The 
solutions with the + sign approaches the regular A = 1 branch in the limit A — > 1, while the solution with the — sign 
is not analytic in that limit. 

As possible continuations of this work, it may be interesting to study solutions where the detailed balance condition 
is softly broken, or where the action is extended to the so called "healthy version" of Hofava gravity. Moreover, the 
study of stability under perturbations may be interesting, as well as the inclusion of a Lifshitz scalar field with a 
symmetry breaking potential that could provide a topologically conserved charge. 
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